A Josephson Junction Microscope for Low-frequency Fluctuators 
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The high-Q harmonic oscillator mode of a Josephson junction can be used as a novel probe 
of spurious two-level systems (TLSs) inside the amorphous oxide tunnel barriers of the junction. 
In particular, we show that spectroscopic transmission measurements of the junction resonator 
mode can reveal how the coupling magnitude between the junction and the TLSs varies with an 
external magnetic field applied in the plane of the tunnel barrier. The proposed experiments offer 
the possibility of clearly resolving the underlying coupling mechanism for these spurious TLSs, an 
important decoherence source limiting the quality of superconducting quantum devices. 



Superconducting quantum circuits have been inten- 
sively tested in various regimes in the past few years, 
from superconducting qubits demonstrating long coher- 
ence times, to superconducting transmission line cav- 
ities coherently coupled to a Single Cooper Pair box 
Such circuits are extremely sensitive to 
very small quanta and defect states, and hence have the 
ability to detect individual microwave photons, charged 
quasiparticles, as well as spurious TLSs within or near 
Josephson junction tunnel barriers @, @, @, Ho, L!Jj. In 
recent experiments 1(1 U|, TLSs were identified through 
spectroscopic measurements of a superconducting phase 
qubit appearing as 'gaps" or "splittings" in the energy 
spectrum. 

The TLS defects can be an unwanted source of de- 
coherence for superconducting quantum bits. The low- 
frequency noise, which has been shown to be a seri- 
ous source of decoherence for superconducting qubits 
[Tol Ull [l2l H3L [3], is very probably induced by such 
amorphous fluctuators inside or near Josephson junc- 
Understanding the origin of these spurious 
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tions 

TLSs, their coherent quantum behavior, and their con- 
nection to ubiquitous 1 /f noise is hence a challenge that 
will be crucial to the future of superconducting quantum 
devices. The behavior of a distribution of these TLSs 
were studied theoretically in 
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Recently, it 

was proposed that TLSs can be viewed as qubits them- 
selves [20j, given their relatively long coherence times. 
However, the microscopic origin and the coupling mech- 
anism between the TLSs and the junction remains un- 
resolved. Generally considered to be connected to the 
amorphous nature of the tunnel barrier [2l|, movement 
of unrestrained atoms or charges may lead to a number 
of possible coupling mechanisms. As originally proposed 
in Ref. [ll|, fluctuations of the TLS could lead to varia- 



tions of the junction critical current. Another possibility, 
requires that the TLSs have fluctuating dipole moments 
which couple to the electric field found within the junc- 
tion tunnel barrier [lQj |. 

Here, we present a scheme that can resolve a variety of 
properties of the TLSs and distinguish between these two 



suggested coupling mechanisms through the use of an ap- 
plied magnetic field. Consider a Josephson junction res- 
onator operating as a high-Q, nonlinear cavity mode [22I ]. 
coupling to a TLS through its canonical phase (or mo- 
mentum) operator. This forms a cavity QED system with 
the junction resonator mode as the cavity and the TLS 
as the atom [23|, HiJ . The junction resonator acts as a mi- 
croscope for studying the behavior of the TLS. We will 
show that microwave transmission in the junction res- 
onator carries spectroscopic, interaction (coupling), and 
spatial information of the TLS. In particular, the junc- 
tion resonator can resolve the coupling mechanism be- 
tween the TLS and the junction. When the TLS couples 
to the junction through the junction critical current, the 
magnitude of the coupling will be strongly modulated 
by changing the strength of the magnetic field oriented 
along the plane of the tunnel barrier. On the contrary, if 
the TLS can only couple through the junction's electric 
field, magnetic fields will have no effect on the magni- 
tude of this coupling. Changes in the coupling magni- 
tude can be observed by measuring the microwave trans- 
mission through the junction resonator. To demonstrate 
this quantitatively, we calculate the transmission and its 
noise spectrum in the "bad cavity" limit [25[, where the 
dissipation of the junction resonator is much faster than 
that of the TLS. Our calculations show that the reso- 
nances in transmission and the noise spectrum strongly 
depend on the coupling strength for chosen TLS. Mean- 
while, the energy distribution, dynamic, dissipative, as 
well as spatial properties of the TLSs can be obtained 
from the measured transmission as well. 



Magnetic field modulation of the coupling. The circuit 
is shown in Fig. [IJ where a Josephson junction with the 
energy Eji together with a superconducting loop forms 
an RF SQUID enclosing a bias flux $5. Magnetic field 
applied along the plane of the tunnel barrier, inside the 
junction [ifj, results in a flux of $1 = At $5 = 

0, the effective Josephson energy of the circuit can be 
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tunnel barrier 



FIG 

<J>i = -^gifi and loop flux: <Ev The Josephson junction has 
energy Ej and capacitance Cj. 



1: Circuit. RF-SQUID loop with in-plane magnetic flux 

h 



that the applied magnetic field $1 shifts the frequency 
of the junction resonator mode which can be tuned by 
controlling the loop flux $5 (see below). 

The driving on the junction resonator can be obtained 
by capacitively coupling the junction to a microwave 
source of 2e cos cJdt(a + a)) with frequency u d and am- 
plitude e. In the rotating frame, the total Hamiltonian 
for the coupled TLS and junction resonator system is 

H t = A c a t a+^cr 2 +^((/9i)(cr + a + cr_a t )+e(a + a t ) (3) 



written as — Ej cos(£ + <pi/2) with 



Ej = Eji 



sin y?i/2 
¥>i/2 



modified by the magnetic field and S being the phase 
difference across the junction. Given a total capacitance 
Cj, the junction behaves as an harmonic oscillator for the 
phase variable 5 with a frequency uo c = yAe^Ej/WCj. 
Consider a TLS inside the junction with an energy hw a 
close to hw c . When the TLS couples with the critical 
current of the junction, the coupling can be derived as 



Jo 



dx cos(£ + <pi — )j d • af(x 



where L is the length of the junction along the direction 
perpendicular to $1 in Fig. HJ 3d describes the polariza- 
tion and magnitude of the coupling, S = (cr x , a y , a z ) are 
the Pauli matrices of the TLS, and f{x — rj) describes the 
spatial profile of the TLS centered at rd [H|,[l3]. For sim- 
plicity, we assume f(x-r d ) = S(x-r d ) and j d = (j x , 0, 0). 
The coupling becomes —Ej\j x cos(5-\-pir d /L)a x . To the 
lowest order of the shifted phase variable (5 + ^i/2), the 
coupling can be written as H c = gdO~ x (d + af) with 



Qd = Ejij x A 



1, 



/ 2e 2 . r d 
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and a and a) being the annihilation and creation opera- 
tors of the phase variable. The coupling hence oscillates 
with and is strongly affected by the applied flux 3>]_. 

In contrast, when the TLS couples to the dielectric field 
within the Josephson junction, the coupling would be 
— 2 q ^f, where p$ is the conjugate of the phase variable 
5, do is the size of the dielectric dipole, and ho is the 
thickness of the junction tunnel barrier. The resulting 
coupling can be written as H c = —ig c cr x (cL — ftf) with 



^0 
h 



e 2 huj c 
2C 



(2) 



not depending on the applied magnetic flux $1. There- 
fore, probing the dependence of the coupling on the mag- 
netic flux will clearly determine which physical mecha- 
nism is more visible even if both effects are present. Note 



with A c = uj c - uj d , A a = uj a - uj d , and g(tpi) given 
by either Eq. ([T]) or Eq. ([2]) depending on the cou- 
pling mechanism. The environmental noise can play 
an important role in the stationary state of the cou- 
pled system. We treat the noise by the Lindblad form 
nC(a)p + 7^£((T_)p + j p C(cr z )p in the master equation, 
which includes dissipation of the junction resonator with 
the rate ft, decay of the TLS with the rate 7^, and de- 
phasing of the TLS with the rate j p . 

The model discussed above describes a typical cavity 



QED system |23|, |24[. The junction resonator mode acts 
as a high-Q cavity driven by microwave source and the 
TLS acts as a two-level atom coupled to the resonator 
mode with a JC-type of interaction. Because of the cou- 
pling, the transmission through the junction resonator is 
imprinted by the properties of the TLS. As we will show 
below, measurement of the microwave transmission pro- 
vides an effective probe, or a microscope, for the TLS. 
Note that the frequencies of the TLSs were observed to 
be separated by 200 MHz on average [10]. With a cou- 
pling magnitude of g ~ 10 MHz, it is reasonable to as- 
sume that only one TLS satisfies the near-resonance con- 
dition: \uo a — uo c \ ~ g and affects the transmission in the 
junction resonator significantly. The TLSs that are far 
off-resonance from the junction resonator mode induce 



(1) small ac-Stark shifts on the order of g /\u a 



Microwave transmission in the junction resonator. To 
quantitatively illustrate the effect of the TLS coupling on 
the transmission in the junction resonator, we study the 
above system in the "bad cavity" limit with 7d,7p <C ft, 
i.e. the dissipation of the junction resonator is much 
faster than that of the TLS [25[. As a result, the junction 
resonator mode adiabatically follows the dynamics of the 
TLS which is governed by the Bloch equation: 



d(a) 
dt 



(4) 



with the Pauli matrices a = (a z , cr + , a_) T . The dynamic 
matrix is 



-71 

in: 



-in* 

iA — 72 






-iA - 72 
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FIG. 2: Transmission in the junction resonator with the cou- 
pling strength ga = (dotted curve), ga = 10 (dashed curve), 
and ga = 20 (solid curve) respectively versus the driving fre- 
quency cod- The parameters are oo c = 7 GHz, A ac = 30 MHz, 
<y d = 0.2 MHz, 7 P = 0.1 MHz, k = 5 MHz, and e = 1 MHz for 
a TLS at the position rd — 0.2. 



and the offset vector is B = A 2 1 • (—71, 0, 0) T with 
i2ge 



k — iA c 



A = A 



a ft 2 + A 2 

Id . , # 2 ft 

72 = y+7p- 



ft 2 + A? 



(5) 
(6) 

(7) 



and 71 = 2(72 — 7 P ) as the dressed parameters in the 
Bloch equation. The coupling between the TLS and the 
junction resonator mode modifies the detuning A and 
the decoherence 71,2 of the TLS, and results in an ef- 
fective driving Q r on the TLS that depends linearly on 
the driving magnitude. The stationary state of the TLS 
is (a) ss = —B. Defining n = Q r /(A + 272), we have 
B z « 1 — 1 77 1 2 ^2 / Ti an d B+ = B*_ w 77/2 under the weak 
driving condition \r]\ <C 1. 

Below, we derive the transmission t{uOd) = \n{a) ss /e\ 2 
at the driving frequency uod, depending on the ratio of 
the output field (a) ss = (e + g(cr_) ss )/(ft + iA c ) and the 
driving field. In Fig. O we plot the transmission versus 
the driving frequency at various values of applied mag- 
netic flux for the coupling magnitude g((fi) = ga defined 
in Eq. (pQ). Let uo^ 2 be the position of the resonance 
peaks in the transmission. At g& = with no coupling, 
a resonance peak appears at uod = 00 c . At (fi = 7r/2, the 
coupling increases to g& = 10 MHz and two resonance 
peaks appear with the separation uo 2 — uo^ = 38 MHz. At 
(fi = 7r, the coupling is ga = 20 MHz and the separation 



of the two resonance peaks becomes uo 2 



50 MHz. 



To explain this result, we make the following approxima- 
tion, 



t(u>d) 



(A a A c -^)2 + „2 A 2- 



(8) 



under the condition 7^, j p <C gi, ft, |A a |, |A C |. With 
no coupling at gd = 0, t(oOd) — ft 2 /(ft 2 + A 2 ) with one 
resonance peak at uod = w c . With finite coupling, the 
approximated transmission in Eq. (|HJ) has two resonance 
peaks with frequencies satisfying A a A c — g 2 = 0., i.e. 
ujp 1 ' 2 -' = (uj a + u c ± \/^ac + 4g 2 )/2 respectively, and 
A ac = ou a — uo c . This shows very good agreement with 
the plots in Fig. O For the coupling given by Eq.([2]), 
the coupling and hence the transmission are not affected 
by the applied magnetic flux. Note that in Eq. ([8]) the 
transmission reaches a minimum at uod — with t(co a ) 
approaches zero. The position of this minimum can be 
used to determine the frequency of the TLS. 

Noise spectrum of the junction transmission. The time 
correlation function of the junction transmission can be 
measured to study the noise spectrum of the transmis- 
sion, and subsequently the properties of the TLS. Let 
X(t) = (a + aJ)(i) be the linearized operator of the 
phase variable (5 + </?i/2). The symmetrized noise spec- 
trum for the operator X(t) can be written as Fxx{u>) — 
+ - In the "bad 

cavity" limit [25[, we have, 



d 
dt 



a a (t) = -(ft - iaA c )a a (t) + e + gcr a (t) (9) 



with a G { — ,+}, cl- = a and a + = a), which forms a 
linear transformation between the Pauli matrices of the 
TLS and the operators of the junction resonator. The 
time correlation function (X(t)X) (and (XX(t))) can 
then be derived from the time correlation functions of 
the Pauli matrices of the TLS. 



Following the quantum regression theorem [27(, the 
time correlation functions of the operators of the TLS 
can be derived from Eq. (|4]) and the matrix A^. At weak 
driving |^| <C 1, the matrix A2 can be approximated to 
be a diagonal matrix with the elements {A zz , 
(not shown) under a second order perturbation of n (77*). 
The time correlation functions can be then derived as 

Mt)ap) = e A - 1 ((a a a p ) - B p B a ) + BpB a (10) 

and similarly for ((T a ap(t)). Among such time correlation 
functions, the dominant contribution to Fxx(w) is by 
* 1. 

We consider the normalized noise spectrum: F{uj) = 
Fxx(u)(k 2 + (uj — A c ) 2 )/g 2 . It can be calculated that 



F(lj) 



72 



72 



( w _A)2 +7 2 ( W + A)2+ 7 | 



(11) 



to an accuracy of 0(|^| 2 ), directly associated with the 
parameters of the TLS. The frequencies of the resonance 
peaks in the noise spectrum are given by uj)ip = ±A 
as defined in Eq. (|6]). The width of the resonance peaks 
is given by 72 as defined in Eq. (0). By measuring the 
resonances in the noise spectrum F(u), detailed charac- 
terization of the TLS can be achieved. For example, at 
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FIG. 3: Normalized noise spectrum F (T _ (T+ (u)) at g — 20 MHz 
and | eta | <C 1 plotted for various driving frequencies. For 
parameters, see the caption of Fig. [2] 



Conclusions. We present a cavity QED scheme for 
studying the properties of the TLSs in the tunnel barrier 
of a Josephson junction. The high-Q oscillator mode of 
the junction resonator acts as a microscope for probing 
the spectral, spatial, and coupling properties of the TLSs 
by the measured the microwave transmission in the junc- 
tion. In particular, our study shows that the coupling 
mechanism between the junction and the TLS can be de- 
termined by applying a magnetic field in the plane of the 
junction tunnel barrier. 
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for Promoting Science and Technology of Univ. Tokyo. 



the driving frequency ujd = oj c (A c = 0), the position of 
the resonances peaks is given by Unp = ±A ac , reveal- 
ing the frequency uo a of the TLS. At A c ^ but A c ~ ft, 
the position of the resonance peaks strongly depends on 
the coupling strength g. In Fig. G3 we plot F(u) for 
g = 20 MHz over a range of driving frequencies. When 
reducing the coupling strength to g <C ft, the noise spec- 
trum shows two resonance peaks with frequencies linearly 
increasing with the driving frequency and crossing each 
other at u)d = u a . 

Spectroscopic Measurements. The spectrum of the TLS 
covers a broad range of frequency. In order to clearly de- 
tect the variations in the coupling and the position of 
the TLSs, we need to control the resonant frequency of 
the junction resonator oj c to sweep past the individual 
TLSs for each successive set of transmission measure- 
ments. This is achieved by varying the applied flux $5 in 
the RF-SQUID loop to adjust the total Josephson induc- 
tance in the RF-SQUID loop. The effective frequency of 
the junction resonator can then be derived as 

CJ C = i COS ^*®b/$o) | 1 (12) 

y LjCj LfrCj 

where Lj is the intrinsic inductance of the Josephson 
junction when $5 = and is the self-inductance of 
the RF-SQUID bias loop. 

Once a TLS is chosen and the frequency of the TLS is 
close to the frequency of the junction resonator mode, we 
want to keep the frequency uo c from changing during the 
transmission measurements while continuously varying 
the magnetic field $1. The bias flux § b in the RF-SQUID 
loop again provides crucial control of the frequency uo c . 
Adjusting the flux $1 inside the tunnel barrier, oj c is 
shifted as we discussed previously. But by tuning the 
flux $5 in the RF-SQUID loop, this frequency shift can 
be compensated to stay in a constant energy contour. 
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